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We present the distance priors that we have derived from the 2015 Planck data, and use these 
in combination with the latest observational data from Type la Supernovae (SNe la) and galaxy 
clustering, to explore the systematic uncertainties in dark energy constraints. We use the Joint 
Lightcurve Analysis (JLA) set of 740 SNe la, galaxy clustering measurements of H{z)s and Da{z)/s 
(where s is the sound horizon at the drag epoch) from the Sloan Digital Sky Survey (SDSS) at 
a = 0.35 and z = 0.57 (BOSS DR12). We find that the combined dark energy constraints are 
insensitive to the assumptions made in the galaxy clustering measurements (whether they are for 
BAO only or marginalized over RSD), which indicates that as the analysis of galaxy clustering data 
becomes more accurate and robust, the systematic uncertainties are reduced. On the other hand, 
we find that flux-averaging SNe la at 2 > 0.5 significantly tightens the dark energy constraints, 
and excludes a flat universe with a cosmological constant at 68% confidence level, assuming a dark 
energy equation of state linear in the cosmic scale factor. Flux-averaging has the most significant 
effect when we allow dark energy density function X(z) to be a free functions given by the cubic 
spline of its value at 2 = 0, i, |, 1; the measured X(z) deviates from a cosmological constant at 
more than 95% confidence level for 0.4 < z < 0.7. Since flux-averaging reduces the bias in the SN 
distance measurements, this may be an indication that we have arrived in the era when the SN 
distance measurements are limited by systematic uncertainties. 
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I. INTRODUCTION 

We continue to search for the unknown cause for the 
observed cosmic acceleration 0 , 0 , a.k.a., dark energy.^ 
Current observational data offer tantalizing hints for de¬ 
viations from a cosmological constant in a simplistic com¬ 
bination of all observational data without critical anal¬ 
ysis. In order to arrive at robust constraints on dark 
energy, we must carefully examine all the data sets sep¬ 
arately, and jointly. One complication is the difficulty to 
detect and model unknown systematic uncertainties in 
the data used for the analysis. 

In this paper, we explore the existence of unknown 
systematic uncertainties by critically analyzing the latest 
observational data from Type la Supernovae (SNe la) 
and galaxy clustering, with the help of distance priors 
from cosmic microwave background anisotropy (CMB) 
data. We use the Joint Lightcurve Analysis (JLA) set 
of 740 SNe la, galaxy clustering measurements of H{z)s 
and Da{z)/s (where s is the sound horizon at the drag 
epoch) from the Baryon Oscillation Spectroscopic Survey 
(BOSS) at z = 0.35 and z = 0.57, and the distance priors 
that we have derived from the 2015 Planck data. 


* email: wang@ipac.caltech.edu 
^ For recent reviews, see 


We describe our method in Sec.II, present our results 
in Sec.Ill, and conclude in Sec.IV. 

II. METHOD 

We focus on exploring the unknown systematic un¬ 
certainties in the current SN la and galaxy clustering 
data using geometric constraints only, with distance pri¬ 
ors from the 2015 Planck data to help tighten parame¬ 
ter constraints. For a conservative and transparent ap¬ 
proach, we marginalize over constraints on the growth 
rate of cosmic large scale structure (which are degener¬ 
ate with the geometric constraints [T^, Il3 ) ■ 

We now give the basic formulae that we will use later 
in the paper. The comoving distance to an object at 
redshift z is given by: 

r(z) = |flfc|“^/^sinn[|llfe|^/^r(z)], (1) 

E{z) = Hiz)/Ho 

Jo T/(Z ) 

where sinn(a:) = sin(a:), x, sinh(a:) for ilk < 0, ilk = 0, 
and ilk > 0 respectively. The Hubble parameter, H{z), 
is given by 

77^(2) (2) 

= Hq [Hm(l + z)^ -l- Hj.(l -|- z)^ -f Hfe(l -f z)^ -f HjfV(z)] , 
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where flm + = 1- The dark energy den¬ 

sity function X{z) = px{z)/px{^)- The $7^ term, with 
Q.r = r7m/(l + Zeq) <C {zeq denotes the redshift at 
matter-rediation equality), is usually omitted in dark en¬ 
ergy studies at z <C 1000, since dark energy should only 
be important at late times. For comparison with the 
work of others and to provide a reference for future sur¬ 
veys, we consider a dark energy equation of state linear 
in the cosmic scale factor a [14) : 

wxia) = Wo + {1 - a)wa. (3) 

In addition, we consider an alternative two parameter 
parametrization of wxia), using wq and wq.o = wx{z = 
0.5): 



= 3^0.5 - 2uio + 3(uio - wo.o)a (4) 

Note that a = | for z = 0.5. It has been shown that 
(wqjWo.s) are significantly less correlated compared to 
(wo,Wo) [il]- 

Finally, we consider a model-independent parametriza¬ 
tion of X(z), where X(z) is a free function of redshift 
given by the cubic spline of its value at z = 0, |, |, and 
1. We assume that X{z > 1) = X{z = I). 


A. CMB data 


We use CMB data in the condensed form of the CMB 
shift parameters (a.k.a., distance priors) M- 


R = ^JnmH^r{z^)/c, 

la = 7rr(z*)/rs(z,). (5) 


These two parameters, R and la, together with Wf, = 
ilbh^, provide an efficient summary of CMB data as far 
as dark energy constraints go [H, jl^ ■ 

CMB data give us the comoving distance to the 
photon-decoupling surface r(z*), and the comoving sound 
horizon at photo-decoupling epoch rs(z*) [l^. The co¬ 
moving sound horizon at redshift z is given by 


Cs dt' 


rs[z = 


= cH, 


-1 


dz' 


Eiz'Y 


= cH, 


-1 


da' 


0 J3{1 + Rba')a''^E^{z') 


( 6 ) 


where a is the cosmic scale factor, a = 1/(1 -I- z), 
and a'^E^{z) = ^m{a + Oeq) + VLkO^ -\- flxX{z)a‘^, 
with Oeq = f7rad/I7m = 1/(1 + ^eq), and Zeq = 

2.5 X 10'^nm^^(Tc'MB/2.7K)“^. The sound speed is 

Cs = 1/+ Rb a), with Rba = ^pb/{ApY, Rb = 
31500I7b/i2(j^^^/2.7K)-4. We take Tcmb = 2.7255. 


The redshift to the photon-decoupling surface, z*, is 
given by the fitting formula [l^ : 

z* = 1048 [1 -k 0.00124(I7b/i2)-°-^3®] [l -f gi{Vtmh^YY , 

(7) 

( 8 ) 

(9) 


where 


-0.238 


5i = 


52 = 


0.0783 (flfc/i^) 

TT39ffi7lW)o^ 

0.560 

l-k21.1(I7b/l2)l-81 


The redshift of the drag epoch Zd is well approximated 

by HI 

[1+0(0,^)“]. (10) 


1 -k 0.659(17^/12)0.828 


where 


hi = 0.313(I7„/i 2)-° ‘‘^® [l-k0.607(I7™/i2)0-674j ^(11) 

62 = 0.238(I7„/i2)‘’'22^ ( 12 ) 

Since the constraints on {la, R,ujb,ns) are not sensi¬ 
tive to the assumption about dark energy (^ . we are 
able to use the Planck archiv to obtain constraints on 
{la,R, 0 Jb,ns) from the 2015 Planck data. We use data 
from the Planck archive that include both temperature 
and polarization data, as well as CMB lensing. As we 
have shown in earlier work [l^ . the one dimensional 
marginalized probability distributions of {la, R,ujb,ns) 
are well fitted by Gaussian distributions. For the Planck 
2015 data, {la, R,Wb,ns) are given by Gaussian distribu¬ 
tions with the following means and standard deviations, 
without assuming a flat Universe: 

{la) = 301.76, cr(/a) = 0.093 

(i?) = 1.7474, cr(A) = 0.0051 

{wb) = 0.02228, cr(wb) = 0.00016 

{ns) = 0.9659, cr(n^) = 0.0048 (13) 

with the normalized covariance matrix of {la, R,Wb, Us): 


( 1.0000 0.4529 -0.3507 -0.3576 \ 

0.4529 1.0000 -0.7000 -0.7780 | 

-0.3507 -0.7000 1.0000 0.5296 

V -0.3576 -0.7780 0.5296 1.0000 / 


Assuming a flat Universe, the Planck 2015 data give 
{la, R,ujb,ns) well fit by Gaussian distributions with the 
following means and standard deviations: 

{la) = 301.77, cr(/a) = 0.090 

(i?) = 1.7482, ct(A) = 0.0048 

{ujb) = 0.02226, a{wb) = 0.00016 

{ns) = 0.9653, cr(n^) = 0.0048 (15) 

with the normalized covariance matrix of {la, R,Wb, ns): 

/ 1.0000 0.3996 -0.3181 -0.3004 \ 

0.3996 1.0000 -0.6891 -0.7677 | , , 

-0.3181 -0.6891 1.0000 0.5152 

V -0.3004 -0.7677 0.5152 1.0000 / 
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We have included Us in our distance priors for com¬ 
pleteness. For the remainder of this paper, we marginal¬ 
ize the CMB distance priors over n^,. This means drop¬ 
ping the 4th row and 4th column from the normalized 
covariance matrix of (/a, .R, Wf,, ris), then obtain the co- 
variance matrix for {la,R,uJb) as follows: 

Co\cMB{Pi,P]) = a{pi) a{pj)'^OYmCoYcMB{Pi,Pj), 

(17) 

where i,j = 1, 2, 3. The rms variance cr(pi) and the nor¬ 
malized covariance matrix NormCovcMB are given by 
Eas. lfT^ and (fldll without assuming a flat universe, and 
Eas. (fT5|) and (fT6|l for a flat universe. 

We include the Planck distance priors by adding the 
following term to the of a given model with pi = la, 
P 2 = i?,and P3 = LOb- 

XCMB = ^P^ [CovcMg{pi,pj)] Apj, Ap, = 

(18) 

where are the mean from Eq. (TT^ (without assum¬ 
ing a flat universe) and Eg. (1X51) (assuming a flat uni¬ 
verse), and CovQ\^g is the inverse of the covariance ma¬ 
trix of [la, R, ujb] from Eg. lflTl) . Note that pi = Ug should 
be added if the constraints on Ug are included in the 
galaxy clustering data. 


B. Analysis of SN la Data 


The distance modulus to a SN la is given by 


po = m — M = 5 log 


dLjz) 

Mpc 


+ 25, 


(19) 


where m and M represent the apparent and absolute 
magnitude of a SN. The luminosity distance (Il^z) = 
(1 + z)r{z), with the comoving distance r{z) given by 
Eq.lP). 

We use the JLA set of 740 SNe la processed by Betoule 
et al. (2014) [2^. They give the apparent B magnitude, 
niB, and the covariance matrix for Am = ms — Wmod, 
with [2^ 


rrimod = 51ogio7?L(^|s) - aXi +/3C + M, (20) 


where 'Db{z\s) is the luminosity distance multiplied by 
Hq for a given set of cosmological parameters {s}, Xi 
is the stretch measure of the SN light curve shape, and 
C is the color measure for the SN. A4 is a nuisance pa¬ 
rameter representing some combination of the absolute 
magnitude of a fiducial SN la, M, and the Hubble con¬ 
stant Hq. Xi is assumed to be different for SNe la with 
different host stellar mass: 


A4 = Ml for host stellar mass < IO^^Mq 
M = M 2 otherwise (21) 

Since the time dilation part of the observed luminosity 
distance depends on the total redshift Zhei (special rela¬ 
tivistic plus cosmological), we have 

'Dl{z\s) = + 2:hei)r(z|s), (22) 


where 2 and Zhei are the CMB restframe and heliocentric 
redshifts of the SN. 

For a set of N SNe with correlated errors, we have 

= Am^ • C~^ ■ Am (23) 

where Am is a vector with N components, and C is the 
N X N covariance matrix of the SNe la. 

Note that Am is equivalent to Apo, since 


Am = mB — rn^nod = [mB + aXi — (5C] — Xi. (24) 

The total covariance matrix is 

C — Dgtat + Clgtat + Cgys, (25) 

with the diagonal part of the statistical uncertainty given 

hyUm 

Ti _L J_ n-2 _L 

-L'stat.ii ' ^int ^lensing 

+a^cr^i j + P'^O'c.i + ‘^CnCmBXi.i ~ ‘^PCrriBC.i 

-2apCx^c,^, (26) 

where CmBXi,i, CmBCy, and CxiCy are the covariances 
between tob, Ai, and C for the i-th SN. Note the Betoule 
et al. (2014) included host galaxy correction in Cgtat + 
Cgys (see Eq.(ll) of [H). 

The statistical and systematic covariance matrices, 
Cgtat and Cgys, are generally not diagonal [1^, and are 
given in the form: 

Cgtat + Cgys = VQ + a^Va+0^Vb + 2aV[)a — 2liV[)b — 2aj3Vab- 

(27) 

where Fq, K, 14, Voa, Vob, and Vab are ma¬ 
trices given by Betoule et al. at the link 

http://supernovae.in2p3.fr/sdss_snls_jla/ReadMe.htmb 
Cgtat includes the uncertainty in the SN model. Cgyg 
includes the uncertainty in the zero point. Note that 
Cgtat and Cgys do not depend on Xi, since the relative 
distance moduli are independent of the value of Xi [2^ . 

We refer the reader to Conley et al. (2011) [2^ and 
Betoule et al. (2014) for detailed discussions of the 
origins of the statistical and systematic errors. 

In order to explore the existence of unknown system¬ 
atic effects, we apply flux-averaging to the JLA SNe la at 
z > 0.5. Flux-averaging was proposed to reduce the sys¬ 
tematic bias in distance measurement due to weak lensing 
magnification of SNe la [2§ - [30| : it has the additional ben¬ 
efit of reducing the bias in distance estimate due to other, 
possibly unknown systematic effects [3l| . This is because 
flux-averaging effectively reduces a global systematic bias 
into a local bias with a much smaller amplitude, which 
in turn results in a reduced impact on global parameter 
constraints. Since weak lensing does not have a signifi¬ 
cant effect on SN la data (see, e.g., [H), any systematic 
biases in the current SN la data are likely dominated by 
other, presently unknown sources. 

Here we apply flux-averaging in the minimal approach 
of flux-averaging the SNe la in each redshift bin at higher 
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z, and then use the usual “magnitude statistics” (instead 
of “flux statistics” in computing since the 

JLA SNe la have measurement and modeling errors that 
have been effectively Gaussianized in magnitudes. 

For statistics using MCMC or a grid of parameters, 
here are the steps in flux-averaging |2l| in application to 
the JLA SNe la: 

(1) Convert the distance modulus of SNe la into 
“fluxes”, 

F{zi) = io-(Mr“Gd-25)/2.5 ^ f ^ (28) 

\ Mpc / 

(2) For a given set of cosmological parameters {s}, ob¬ 
tain “absolute luminosities”, {C{zi)}, by removing the 
redshift dependence of the “fluxes”, i.e., 

C{zi) = di{zi\s) F{zi). (29) 

(3) Flux-average the “absolute luminosities” {£]} in 

each redshift bin i to obtain < £ : 


Ni 




Z=1 


1 

(30) 

* 1=1 


(4) Place £* at the mean redshift Zi of the i-th redshift 
bin, now the binned flux is 

F{z,)=T/dl{z,\s). (31) 

with the corresponding flux-averaged distance modulus: 

-data(-) ^ _2.51ogioF(z,) + 25. (32) 

(5) Compute the covariance matrix of 7l(zi) and Jlizj): 

Cov [-p{zi),Jl{zj)\ (33) 

1 

Ni Nj 

m— 1 

where is the covariance of the 

measured distance moduli of the Z-th SN la in the Ath 
redshift bin, and the m-th SN la in the j-th redshift bin. 
C{z) is defined by Eas. (l28l) and (l2^ . 

(6) For the flux-averaged data, {jl{'Zi)}, compute 

= X] Am(^z) Cov"^ [■p{zi),Jl{zj)] A7Z(zj) (34) 
u' 

where 

Ajl{zi) = Jl{zi) - fjF{zi\s), (35) 

and 

JlP{z^) = -2.5 logio FP{zi) + 25. (36) 

with FP{zi\s) = ((ii(z|s)/Mpc)~^. 

For the sample of SNe we use in this study, we flux- 
averaged the SNe with dz = 0.04. 


C. Galaxy Clustering Data 

For GC data, we use the measurements of Xh{z) = 
H{z)rs{zd)/c and Xd{z) = DA{z)/rs{zd), where H{z) is 
the Hubble parameter, Da{z) is the angular diameter 
distance, and rs{zd) is the sound horizon at the drag 
epoch. It has been shown that Xh{z) and Xd{z) are more 
tightly constrained by data, and less sensitive to mod¬ 
eling assumptions, compared to F[{z) and Da{z) [^ . 
We use the Xh{z) and Xd{z) measurements from the two- 
dimensional power spectrum measured at z=0.32 and 

z=0.57 from BOSS DR12 galaxies [H, Converting 

the results in [H, [2^ to the same definitions used in this 
paper, we find that for BAO only [2^ 

Xft(0.32) = iL(0.32)rs(zd)/c = 0.0397 ± 0.0021 
a;d(0.32) = £iA(0.32)/rs(zd) = 6.49 ±0.16 

rM(0.32) = 0.41 (37) 

Xh{QAl) = H{Q.bl)rs{zd)/c = 0.0498 ±0.0013 
Xd{QA7)=DA(J)A7)/r,{zd) = 9.18 ±0.13 

rM(0.57) = 0.47 (38) 

For BAO measurements marginalized over BSD [2^, we 
find 


Xh{Q.2>2) = H{QA2)rs{zd)lc 
Xd{0.32) = DA{0.32)/rs{zd) 
rhd{0-32) 


0.0391 ±0.0019 
6.185 ±0.185 
0.5 (39) 


Xh{0.57) = H{0.57)rs{zd)/c 
Xdi0.57)=DAi0.57)/rsizd) 
rhdiOA7) 


0.0476 ±0.0015 

9.18 ±0.15 

0.53 (40) 


Galaxy clustering data are included in our analysis by 
adding Xgc = Xgci + XgC 2 > ^GCi = 0.35 and 
ZGC 2 = 0.57, to the x^ of a given model. Note that 

XGCi = [Cgc (Pi ’Pj)] ^Pj > Api =p^- , 

(41) 

where pi = H{zGCi)rs{zd)/c and p 2 = DAizGGi)lrs{zd), 
with i = 1, 2. 


III. RESULTS 

We perform a MCMC likelihood analysis [s^ to obtain 
0(10®) samples for each set of results presented in this 
paper. We assume flat priors for all the parameters, and 
allow ranges of the parameters wide enough such that 
further increasing the allowed ranges has no impact on 
the results. We constrain dark energy and cosmological 
parameters (wq, Wa, ^m, Xlk, h,u!h), where uib = In 

addition, we marginalize over the SN la nuisance param¬ 
eters {a,/3, Ml, M 2 }. 
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FIG. 1: Marginalized probability distributions of parameters from 
JLA SNe, galaxy clustering data at z = 0.35 and z = 0.57, and 
Planck 2015 distance priors. The solid and dotted curves corre¬ 
spond to using H{z) and Da[z) measurements from BAO only 
measurements, and those from RSD marginalized measurements. 


A. Constrains on wo and Wa 

FigH] shows the marginalized probability distributions 
of parameters from JLA SNe, galaxy clustering data at 
0 = 0.35 and z = 0.57 [H,[2^, and Planck 2015 distance 
priors presented in this paper (see Sec.IIA). The solid and 
dotted curves correspond to using JI(z) and Da(,z) mea¬ 
surements from BAO only measurements, and those from 
RSD marginalized measurements. Figl2] shows the joint 
68% and 95% confidence contours for {wa, wq) and {wa, 
Dfc) corresponding to FiglU with the same line types. 
The combined dark energy constraints seem insensitive 
to the assumptions made in the analysis of galaxy clus¬ 
tering data. 

FiglS] shows the impact of flux-averaging SNe la on 
the marginalized probability distributions of parameters 
from the combination of the same data sets as in FiglT] 
The solid and dotted curves correspond to using SNe la 
with and without flux-aver aging. Figl4] shows the joint 
68% and 95% confidence contours for (wq, wq) and {wa, 
Dfc) corresponding to FigEl with the same line types. 
Clearly, flux-averaging significantly tightens the dark en¬ 
ergy constraints. This may be due to the reduction in 
the bias of distance measurements from flux-averaging, 
which increases the concordance of the data, resulting in 
tighter constraints. 


B. Constraints on wo and wo.s 

Figsl5][8] are similar to FigslTJdl but for parametrizing 
the linear dark energy density uisng wq and W 0.5 (see 
Eq.U), instead of the usual wq and Wa- Fig (S] shows 
the impact of the galaxy clustering analysis technique 
on the marginalized probability distributions of param¬ 
eters from JLA SNe, galaxy clustering data at z = 0.35 
and z = 0.57 [ 2 ^ [ 2 ^, and Planck 2015 distance pri¬ 
ors presented in this paper (see Sec.IIA). The solid and 




FIG. 2: The joint 68% and 95% confidence contours for (wa, 'Wo) 
and (iCa, fifc) corresponding to Fig^ with the same line types. 






0.022 0.023 

O^h* 


FIG. 3: Marginalized probability distributions of parameters from 
JLA SNe, galaxy clustering data at 2 : = 0.35 and 2 : = 0.57 (BAO 
only), and Planck 2015 distance priors. The solid and dotted curves 
correspond to using SNe la with and without flux-averaging. 


dotted curves correspond to using i7(z) and Da{z) mea¬ 
surements from BAO only measurements, and those from 
RSD marginalized measurements. Figl6] shows the joint 
68% and 95% confidence contours for (wo. 5 ) wq) and 
(wo. 5 , fik) corresponding to FiglU with the same line 
types. Again, the assumptions made in the analysis of 
galaxy clustering data have little impact on the combined 



FIG. 4: The joint 68% and 95% confidence contours for (wa, wg) 
and {wa, life) corresponding to Fig[3] The solid and dotted curves 
correspond to using SNe la with and without flux-averaging. 
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JLA SNe + BOSS DR12 + Planck 2015: BAO only (solid); RSD marginalized (dotted) 



0.65 0.7 0.75 0.022 0.023 

h 



FIG. 5: Marginalized probability distributions of parameters from 
JLA SNe, galaxy clustering data at 2 = 0.35 and 2 = 0.57, and 
Planck 2015 distance priors. The solid and dotted curves corre¬ 
spond to using H(z) and D^( 2 ) measurements from BAO only 
measurements, and those from RSD marginalized measurements. 




FIG. 6: The joint 68% and 95% confidence contours for ’^o) 

and (u;o.5? corresponding to FigO The solid and dotted curves 
correspond to using H{z) and D^( 2 ) measurements from BAO only 
measurements, and those from RSD marginalized measurements. 


FIG. 7: Marginalized probability distributions of parameters from 
JLA SNe, galaxy clustering data at 2 = 0.35 and 2 = 0.57 (BAO 
only), and Planck 2015 distance priors. The solid and dotted curves 
correspond to using SNe la with and without flux-averaging. 



FIG. 8: The joint 68% and 95% confidence contours for (rtio.S) 
and ^k) corresponding to Fig[7| The solid and dotted curves 

correspond to using SNe la with and without flux-averaging. 


dark energy constraints. 

FiglT] shows the impact of flux-averaging SNe la on 
the marginalized probability distributions of parameters 
from the combination of the same data sets as in Fig IS] 
The solid and dotted curves correspond to using SNe la 
with and without flux-aver aging. FigjSj shows the joint 
68% and 95% confidence contours for (ico.s, wq) and 
(wo.s, ^k) corresponding to Fig]?] with the same line 
types. Flux-averaging of SNe makes an even more dra¬ 
matic difference in the joint constraints on (icq, wq.s), 
compared to that of {wo,Wa)-- 

C. Constraints on dark energy density fnnction 

We now examine the dark energy constraints when 
we allow the dark energy density function, X(z) = 
px{z)/px{z = 0), to be a free function, given by the 
cubic spline of its value atz = 0,l/3,2/3,l, and assum¬ 
ing that X{z > 1) = X(z = 1). 

FigEI shows the impact of the galaxy clustering anal¬ 
ysis technique on the marginalized probability distri¬ 


butions of parameters from JLA SNe, galaxy cluster¬ 
ing data at z = 0.35 and z = 0.57 [1^ [l^, and 
Planck 2015 distance priors presented in this paper (see 
Sec.IIA). The solid and dotted curves correspond to using 
H{z) and Da{z) measurements from BAO only measure¬ 
ments, and those from RSD marginalized measurements. 
FigfTOl shows the impact of flux-averaging the SNe la on 
the marginalized probability distributions of parameters 
from the same combination of data sets. The solid and 
dotted curves correspond to using SNe la with and with¬ 
out flux-aver aging. Again, we find that the assumptions 
made in the galaxy clustering data analysis have little 
impact on the combined dark energy constraints, while 
flux-averaging of SNe la has a significant impact on these 
constraints. 

Fig llll shows the dark energy density function X(z) = 
px{z)/pxi0) measured from JLA SNe, galaxy clustering 
data at z = 0.35 and z = 0.57 [H, and Planck 2015 
distance priors presented in this paper (see Sec.IIA). The 
shaded regions indicate the 68% confidence region, while 
the outer envelope indicates the 95% confidence level. 
The densely shaded and sparsely shaded regions corre¬ 
spond to using SNe la with and without flux-averaging. 
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JLA SNe + BOSS DR12 Planck 2015: BAO only (solid); RSD marginalized (dotted) 





FIG. 9: Marginalized probability distributions of parameters from 
JLA SNe, galaxy clustering data at 2 = 0.35 and 2 = 0.57, and 
Planck 2015 distance priors. The solid and dotted curves corre¬ 
spond to using H{z) and D^( 2 ) measurements from BAO only 
measurements, and those from RSD marginalized measurements. 


FIG. 11: The dark energy density function X(z) = px{z)/px{0) 
measured from JLA SNe, galaxy clustering data at 2 = 0.35 and 
2 = 0.57 [^ . l24 | . and Planck 2015 distance priors presented in this 
paper (see Sec.IIA). The shaded regions indicate the 68% confi¬ 
dence region, while the outer envelope indicates the 95% confidence 
level. The densely shaded and sparsely shaded regions correspond 
to using SNe la with and without fiux-averaging, respectively. 



FIG. 10: Marginalized probability distributions of parameters 
from JLA SNe, galaxy clustering data at 2 = 0.35 and 2 = 0.57 
(BAO only), and Planck 2015 distance priors. The solid and dotted 
curves correspond to using SNe la with and without flux-averaging. 


respectively. Flux-averaging has the most significant ef¬ 
fect here — the measured X{z) deviates from X{z) = 1 
{w = —1) at more than 95% confidence level for 0.4 < 
0 < 0.7. 

IV. DISCUSSION AND SUMMARY 

We have explored the existence of unknown systematic 
uncertainties in the current SN la and galaxy clustering 
data, with the help of the latest CMB distance priors. 
We use the JLA set of 740 SNe la from Betoule et a. 
(2014) [ 2 ^, and the measurements of H{z) and Da{z) at 
z = 0.35 and z = 0.57 from BOSS DR12 data by Gil- 
Marin et al. 2016 [H, 1^ . We have derived the CMB 
distance priors from Planck 2015 data, in the form of the 
mean values and covariance matrix of {/a, i?, n-s}, 

which give an efficient summary of Planck data in the 
context of dark energy constraints (see Eos. (11311161) 1. 

It is remarkable that the Planck distance priors that 
we have derived from the 2015 Planck data have un¬ 
certainties that are within 10% of the forecasted errors 
for Planck by Mukherjee et al. (2008) [s^. This indi¬ 
cates that Planck has achieved its forecasted precision in 
cosmological constraints. We note that Huang, Wang, 
& Wang [ 3 ^ independently derived similar but slightly 
different constraints from Planck 2015 data at approxi¬ 
mately the same time. 

We find that the combined dark energy constraints are 



















































































insensitive to the assumptions made in the galaxy cluster¬ 
ing measurements (whether they are for BAO only [2^ or 
marginalized over RSD [13 )i independent of the dark en¬ 
ergy parametrization used (see FiglU Figl^J FigjSl FiglHl 
Fig®. We note that the published BAO only constraints 
in [23| differ from those in the earlier arXiv version, and 
are closer to the RSD marginalized constraints in [^ . 
This is reassuring, as it indicates that as the analysis of 
galaxy clustering data becomes more accurate and ro¬ 
bust, the systematic uncertainties are reduced. 

On the other hand, we find that flux-aver aging SNe 
la at z > 0.5 significantly tightens the dark energy con¬ 
straints, and excludes w = — 1 at greater than 68% confi¬ 
dence level (see Fig® Fig® Fig® Fig® FigdOl FigHT]). 
Flux-averaging has the most significant effect when we 
allow dark energy density function X{z) = px{z)/px{^) 
to be a free functions given by the cubic spline of its 
value at z = 0, i, |, 1; the measured X{z) deviates from 
A(z) = 1 (w = —1) at more than 95% confidence level 
for 0.4 < z < 0.7 (see Fig lTTl) . This is somewhat surpris¬ 
ing, since for SN data with redshift-dependent systematic 


biases that are negligible compared to statistical errors, 
flux-averaging of SNe should give somewhat less stringent 
constraints on dark energy |28| . Since flux-averaging re¬ 
duces the bias in the SN distance measurements , this 
may be an indication that we have arrived in the era when 
the SN distance measurements are limited by systematic 
uncertainties. 

Identifying and correctly modeling systematic effects 
will be key in illuminating the nature of dark energy. 
Future dark energy surveys from space [37l - l40l| will be 
designed to minimize systematic uncertainties. We can 
expect dramatic progress in the next decade in our quest 
to shed light on dark energy. 
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